Weak coupling d-wave BCS superconductivity and unpaired electrons in overdoped La 2 - A Sr x Cu0 4 

single crystals 
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The low-temperature specific heat (SH) of overdoped La2_ v Sr r Cu04 single crystals (0.178 < x < 0.290) has 
been measured. For the superconducting samples (0.178 < x < 0.238), the derived gap values (without any 
adjusting parameters) approach closely onto the theoretical prediction Ao = 2A4k B T c for the weak-coupling 
rf-wave BCS superconductivity. In addition, the residual term y(0) of SH at H = increases with x dramati- 
cally when beyond x ~ 0.22, and finally evolves into the value of a complete normal metallic state at higher 
doping levels, indicating growing amount of unpaired electrons. We argue that this large y(0) cannot be simply 
attributed to the pair breaking induced by the impurity scattering, instead the phase separation is possible. 

PACS numbers: 74.25.Bt,74.20.Rp, 74.25.Dw, 74.72.Dn 



I. INTRODUCTION 

For hole-doped cuprates, it is now generally perceived that 
the superconducting state has robust c/- wave symmetry. 1 In the 
underdoped region, due to the presence of the pseudogap and 
other possible competing orders, 2,3 the measured quasipartcle 
gap may not reflect the real superconducting gap. In contrast, 
in the overdoped region, the normal state properties can be de- 
scribed reasonably well by the Fermi liquid picture, 4 although 
still with electronic correlation to some extent. 5 Under this cir- 
cumstance, one may think that the overdoped cuprate provides 
a clean gateway to the intrinsic high-r c superconducting state. 
To accumulate experimentally accessible parameters, such as 
the superconducting gap, and compare them with the mean 
field BCS prediction in this very region is thus expected to be 
particularly valuable. 

Another puzzling point in the overdoped cuprates is that 
the superfluid density p s determined by muon spin relaxation 
(jt/SR) technique decreases just as the transition temperature 
T c when beyond a critical doping point p c ~ 0.19. 6,7,8 This is 
actually not demanded by the BCS theory. The decrease of 
p s , first reported in Tl2Ba2Cu06+,5 (T12201) and subsequently 
confirmed in other families of cuprates, 8,9 was attributed to 
the unpaired carriers at T — > in overdoped cuprates. 6 Sim- 
ilar conclusion was also drawn from studies of the optical 
conductivity 10,11 and magnetization. 12 Recently, the Meissner 
volume fraction was revealed to decrease as T c with increas- 
ing doping in overdoped La2- A Sr A Cu04 (LSCO) and the result 
was suggested to be consistent with the occurrence of a phase 
separation into superconducting and normal-state regions^ It 
is thus highly desired to use the specific heat (SH) which is 
very sensitive to the quasiparticle density of states (DOS) at 
the Fermi level to directly probe these unpaired charge carri- 
ers. 

In this paper we shall address these two issues by the low- 
temperature SH which has established its importance to iden- 
tify the pairing symmetry in high-r e cuprates over the past 



decade. 14 Recently, the quantitative analysis shows that it pro- 
vides a bulk way to obtain the magnitude of the superconduct- 
ing gapji^ 6 - By analyzing the field-induced SH, it is found 
that the pairing symmetry in the overdoped regime (up to 
x = 0.238) is still li-wave and the derived gap values Ao 
approach closely onto the theoretical prediction of the weak- 
coupling (i-wave BCS superconductivity. Our data also reveal 
a quick growing of the residual term y(0) of SH at T — > 
with increasing doping, which cannot be simply attributed to 
the pair breaking induced by the impurity scattering. 



II. EXPERIMENT 



Single crystals of LSCO were grown by the traveling- 
solvent floating-zone method. Details of the sample prepa- 
ration have been given elsewhere ji^Jl The Sr content of the 
sample, x, taken as the hole concentration p, was determined 
from the inductively coupled plasma measurement. Figure 
1(a) shows the dc magnetization curves measured in 10 Oe 
after the zero-field cooling (ZFC) mode using a SQUID mag- 
netometer, where the onset of the diamagnetic signal was de- 
fined as the critical temperature T c . Six samples have been 
measured, with x = 0.178, 0.202, 0.218, 0.238, 0.259, 0.290 
and T c = 36, 30.5, 25, 19.5, 6.5, and below 1.7 K, respec- 
tively. As shown in Fig. 1(b), the T c can be described well 
by the empirical formula 18 T c /T™ x = 1 - 82.6(jc - 0.16) 2 
with T™ 3 " = 38 K. The SH measurements were performed on 
an Oxford Maglab cryogenic system using the thermal relax- 
ation technique, as described in detail previously. 19 The tem- 
perature was down to 1 .9 K and the magnetic field was applied 
parallel to c-axis in the measurements. 
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FIG. 1: (color online) (a): Temperature dependence of DC magneti- 
zation measured in ZFC mode at 10 Oe. The curves are normalized 
to unity with the value at the lowest temperatures, (b): Doping de- 
pendence of T c (squares). The empirical formula (see text) is plotted 
as the solid line. 



HI. RESULTS AND DISCUSSION 

The raw data of SH for all six samples in various H at 
T < 7 K are shown in Fig. 2. To separate the electronic 
SH from other contributions, the data are fit to C(T,H) = 
yT + C P hT + Csch(7\ H), where C p hT = /3T } is the phonon 
SH. Csdi(T,H) is the two-level Schottky anomaly with the 
form D/T 2 in H = and /jcV/(1 + e x f (x = gfi B H/k B T) 
in H + 0, where /i B is the Boher magneton, g the Lande fac- 
tor, and / the concentration of spin- 1/2 particles. The first 
linear- T term, yT, contains the electronic SH and resides in 
the heart of the present study. As demonstrated by the solid 
lines in Fig. 2(a)-(e), all data sets can be described reason- 
ably well by the above expression. For C p h and the Debye 
temperature ©o the values derived here are in consistent with 
previous reports on the sample at similar doping leveli^S For 
Csch, the yielded / is relatively constant at high fields with an 
averaged value ~ 3 mJmoL'K -1 for different samples. This 
low / reflects the small contribution of Cs c h to the total SH 
and assures the reliable determination of y. 

In zero-field, after removing the Schottky term Cs c h and 
by doing a linear extrapolation to the data shown in Fig. 
2(f) to T - K we can determine the residual term y(0) 
of SH. By increasing H, an increase in y is observed for 
0.178 < x < 0.238, as shown in Fig. 2(a)-(d), correspond- 
ing to y = y(0) + y(H) with y(H) the coefficient of the 
field-induced SH. For li-wave superconductors, it was theo- 
retically pointed out that the y(H) is proportional to V77 due 
to line nodes of the gap, 21 which has been confirmed in several 



experiments. 19 Figure 3 summarizes the field dependence of 
the y{H) for the overdoped LSCO. It is clear that for all four 
samples, the y(H) is well described by A V77 with A a doping- 
dependent constant, as exemplified by the solid lines in Fig. 3. 
This indicates that in overdoped LSCO the gap remains robust 
J-wave symmetry. 

Next one can further obtain the gap magnitude by investi- 
gating y{H) quantitatively. Fundamentally, y{H) arises from 
the Doppler shift of the quasiparticle spectrum near the nodes 
due to the supercurrent flowing around vortices and thus di- 
rectly relates to the slope of the gap at the node, va = 2Ao /ftkf 
with Ao the of-wave maximum gap in the gap function A = 
Ao cos(2</>), k F the Fermi vector near nodes (taking ~ 0.7 A -1 
as obtained from ARPES 22 ). Explicitly, the relation between 
Va and the prefactor A is given by 

A - l^t [~*~ nV >™i a fix 
3ft V^o d v A 

where Oo is the flux quantum, n the number of CuCh planes 
per unit cell, d the c-axis lattice constant, V mo i the volume per 
mole, and a - 0.465 for a triangular vortex lattice! 2 ^! Inset 
of Fig. 3 shows the doping dependence of A by fitting the 
data to y(H) = A V77. Thus, with the known parameters for 
LSCO (n = 2,d= 13.28 A, V mol = 58 cm 3 ), the doping de- 
pendence of v A and Ao can be derived without any adjusting 
parameters according to Eq. (1). In this way we extracted the 
gap values A = 9.2 + 0.7, 6.6 ± 0.3 and 5.6 ± 0.3 meV for 
x = 0.178, 0.202 and 0.218, respectively (for x = 0.238, the 
observed A should be corrected due to the volume correction 
which will be addressed later). It can be seen immediately that 
Ao decreases with increasing doping in the overdoped LSCO, 
concomitant with the decrease of T c . The same trend of T c 
and A () with overdoping implies that the suppression of super- 
conductivity mainly comes from the decrease in the pairing 
gap. Figure 4 plots the doping dependence of An, together 
with the value extracted from SH in underdoped and optimal 
doped LSCO single crystals.^ 2 ^ For comparison, the weak- 
coupling li-wave BCS gap relation Ao = lAAksTc is also plot- 
ted as a dotted curve in Fig. 4(a) and a dotted horizontal line 
in Fig. 4(b),™ where T c is determined by the empirical for- 
mula described before. Remarkably, beyond x ~ 0.19, the 
experimental data approaches closely onto the theoretical pre- 
diction, revealing a strong evidence for the weak coupling d- 
wave BCS superconductivity. Previous results about Ao deter- 
mined by scanning tunnelling spectroscopy 27 and penetration 
depth measurements 28 are in excellent quantitative agreement 
with our present results, which strongly supports the validity 
of the present analysis. 

Now we examine the implication of the residual term y(0) 
in zero-field. Figure 5(a) summarizes the doping depen- 
dence of y(0), where the values from previous studies are 
also includedfi^i 2 ^ For comparison, the normal-state SH co- 
efficient yw in the corresponding doping region is shown 
together. 29 We can see that the y(0) increases with doping up 
to x = 0.259. For the nonsuperconducting x = 0.290 sample, 
the y(0) is actually the y#, which shows good consistency with 
the previous value. Note that for x = 0.259, the y(0) is already 
close to the reported y#. Close to the optimal doping point 
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FIG. 2: (color online) Temperature and magnetic field dependence of SH in C/T vs T 2 plot, (a)-(e): Raw data for all six samples (symbols). 
ji H varies up to 12 T for 0. 178 < x < 0.238 while up to 2 T for x = 0.259 and 0.290. The solid lines represent the the theoretical fit (see text). 
The fits are limited to T = 1, 6, 5, and 4 K for x = 0.178, 0.202, 0.218, and 0.238, respectively. For x = 0.259 and 0.290, the fits are ranged to 
7 K. (f): Replot the data at ^i H = T for all samples (symbol: □ = 0.178, O = 0.202, A= 0.218, v = 0.238, = 0.259, <= 0.290). The dotted 
lines are extrapolations of the data down to T = K with the Schottky anomaly subtracted. 



the small y(0) may be attributed to the impurity scattering by 
which a finite DOS is generated for a li-wave superconductor. 
However the large y(0) appearing beyond x ~ 0.22 cannot 
be simply attributed to this reason. This can be understood 
by having an estimation on the impurity scattering induced 
DOS y™ p in the superconducting state, which has the rela- 
tion j^Ijn = (2yo/7rA ) ln(A /yo) with y the pair breaking 
parameter. 30 Also in the unitarity limit, yo - 0.6 VTAo with 
F = 1 /2to the normal state quasiparticle scattering rate which 
can be estimated from the residual resistivity po = m* /ne 2 To 
and the plasma frequency oj p - -^ne 2 /e '«*. With po = 26 
juflcm and a> p ^ 1.2 eV for x = 0.238, 31 one gets ftF ^ 2.5 

meV. Assuming Ao = 2.l4ksT c , we obtain y™ p /yv - 0.2 and 
therefore y™ p ^ 2.9 mJmoF'Kr 2 for x = 0.238, which is far 
below the y(0). 

Furthermore, if the large y(0) is completely induced by the 
impurity scattering, the field-induced y{H) at low H is ex- 
pected to deviate from the V77 dependence and instead show 
an H In //behavior: y(H) a A(H/H c2 )]n[B(H c2 /H)l where 



A = Aoa 2 y;v/8yo with B = n/2a 2 £- In Fig. 6 we present the 
fits to y(H) for H < 4 T with this function. First, we leave A 
and H C 2 as both free fitting parameters (fitl) and the best fit is 
shown in Fig. 6(a). As shown in Fig. 6(b), this yields the pa- 
rameter fx H c2 < 4 T for all samples 0.178 < x < 0.238. The 
rather low H c2 is physically unacceptable. At the same time, 
from the parameter A, the coefficient of the residual specific 
heat y™ p can be calculated using the expressions and the y^ 
described above. It can be seen, for x > 0.218, the obtained 
y|™ p is also inconsistent with the experiment. Secondly, we 
try to fit the data with H c2 fixed within the values shown in 
the shaded region in Fig. 6(e) (fit2). The transport and Nernst 
effect measurements have indicated that fioH c2 ~ l.5T c (H c2 
in Tesla and T c in Kelvin) for the overdoped LSCOi 3 ^ 3 ^ The 
current SH suggests hqH c2 > 12 T for all samples. Hence, 
in Fig. 6(e) the lower limit of the shaded region is set to be 
HoH c2 = 12 T and the upper limit to be fioH c2 = 2T C . In this 
case, we could not obtain a satisfactory fit to the data, as in- 
dicated by the typical result shown in Fig. 6(d). Again, the 



4 



£ 2 

g 
~ 1 



x = 0.178 




FIG. 3: (color online) Coefficient of the field-induced linear-r spe- 
cific heat for 0.178 < x < 0.238, y(H) = y - y(0) at T = K (sym- 
bols). The solid lines are the fits to y(H) = A y[H. Inset: Doping 
dependence of the prefactor A (mJmol~ 1 Kr 2 T~ a5 ). 
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FIG. 4: (color online) Doping dependence of the superconducting 
gap A obtained from SH measurements, (a): A vs x. The dashed 
line is a guide to the eye. (b): Ao/k B T c vs x. The values from Ref. 
16 (half-filled squares) and Ref. 25 (triangle) are also included. The 
weak-coupling d-wave BCS value, Ao = 2.\Ak B T c (Ao/k B T c = 2.14) 
is plotted as a dotted curve (horizontal line) in (a) ((b)). For x > 0.19, 
the experiments are consistent with the BCS prediction. 



FIG. 5: (color online) (a): Doping dependence of the y(0) in zero 
field (filled squares) and the normal state SH coefficient, y N (circles). 
The y(0) from Ref. 19 (half-filled square) and Ref. 25 (up-triangles) 
are also shown. The y N is quoted from Ref. 29. (b): Doping de- 
pendence of the y(0) in (a) normalized by y N , y(0)/yA>. The same for 
Tl-2201 (diamonds) is quoted from Ref. 36. The normalized residual 
spin Knight shift in LSCO, 35 N ics /Nn, is also shown for comparison 
(down-triangles). 



obtained y™ p is contradictory to the measurement (Note, for 
a given sample, one would obtain the lower y™ p with a higher 
H&) (Fig. 6(f)). Therefore, it seems that the impurity scat- 
tering effect could not account for the field dependence of the 
7(H). 

The above analysis suggests that in highly doped LSCO the 
y(0) mainly comes from contributions other than the impurity 
scattering. We attribute it to the presence of nonsuperconduct- 
ing metallic regions. This can be corroborated by simultane- 
ously having a good consistency with the Volovik's relation 
y(H) = A VH and the very large ratio y{Q)ljN on the single 
sample x = 0.238. Figure 5(b) shows the ratio of y(Q)lyN 
together with the normalized residual spin Knight shift, 35 
N ks /Nn, another probe of the residual DOS in the supercon- 
ducting state. In overdoped T12201 the low-temperature SH 
has been measured by Loram et al. 36 and the y(5K)/yAi is plot- 
ted together. We can see that all these quantities show a rapid 
increase with overdoping, indicating a generic property. Actu- 
ally in LSCO previous results also showed the rapid increase 
of y(0) with doping in highly overdoped region although those 
experiments were done on polycrystalline samplesi22i2!i One 
may argue, in LSCO, that there is a high-temperature tetrago- 
nal to orthorhombic structural transition near x a 0.2^i which 
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FIG. 6: (color online) Fit low field y(H) to the HlnH function (see 
text), (a): Fitl with A and H c2 as both free fitting parameters. The 
yielded parameters (A translating to y™ p ) are shown in (b) and (c). 
(d): Fit2 with H c2 fixed within the values shown in the shaded region 
in (e). The dotted and dashed lines in (e) denote H c2 = IT C and l.5T c , 
respectively. The yielded parameter y™ p is shown in (f) as the shaded 
region. For comparison, the experimental y(0) is shown in (c) and (f) 
(circles). 



may induce the rapid increase of y(0). We note that, however, 
in Pr-doped LSCO, 38 this subtle structural transition can be 
tuned to much higher doping level, but the superconducting 
dome remains unchanged, indicating that the hole concentra- 
tion rather than the slight structure distortion plays a domi- 
nant role here. Furthermore, as shown in Fig. 5, a very sim- 
ilar residual y(0) appears in T12201, a system without such a 
structural transition. 

The presence of nonsuperconducting metallic phase implies 
immediately a decrease of the superconducting volume frac- 
tion. This can just explain the field dependence of the SH in 
x = 0.238 and 0.259 samples. For x = 0.238, the observed A 
is even lower than that for x = 0.218, implying a significantly 
reduced superconducting volume fraction. Taking this into 
account, the A used to derive the Ao for this sample, should 
be corrected roughly as Ay^llyN — 7(0)], with the assump- 
tion that the volume ratio is similar to the DOS ratio. The gap 
value yielded with this correction is about 3.5 meV, which also 



scales with T c in if-wave BCS manner and is plotted in Fig. 4. 
For x = 0.259, the sample shows a large y(0) being close 
to the jm, indicating a rather small superconducting volume 
fraction. 

So far, we have shown that in overdoped LSCO, the super- 
conducting gap decreases with increasing x and at high dop- 
ing levels, there exist nonsuperconducting metallic regions at 
T — > 0. Let us discuss the implications of both effects. Pre- 
viously it was suggested that the suppression of T c in over- 
doped regime may come from the increasing pair breaking 
effect. Our present result, however, does not support this pro- 
posal since T c is found to scale with Ao in good agreement 
with c/-wave BCS theory, which implies that the decrease in T c 
should originate from an underlying reduction in the pairing 
strength. This point may also be helpful to elucidate the ori- 
gin of the presence of nonsuperconducting metallic regions in 
the overdoped sample, which is yet unclear. Currently several 
scenarios have been proposed to account for this anomalous 
phenomenon. One is that the overdoped cuprate may spon- 
taneously phase separate into the hole-poor superconducting 
region and the hole-rich normal Fermi liquid region due to 
the competition in energy between these two phases.— An- 
other scenario is associated with the microscopic hole dop- 
ing state. 13 It was speculated that in the overdoped regime the 
holes were doped directly into the Cu3d orbital rather 02p, 40 
which is expected to produce free Cu spins and/or disturb the 
antiferromagnetic correlation between Cu spins. Around these 
holes, the superconductivity is destroyed, forming the normal 
state region. Our present result seems supporting this scenario 
with the assumption that the superconductivity is magnetic in 
origin and the suppression of Ao originates from the disturbing 
of the antiferromagnetic correlation with overdoping. 



IV. SUMMARY 

In summary, low-temperature SH in overdoped LSCO sin- 
gle crystals has revealed two interesting findings: (1) The 
field-induced SH follows the prediction of of-wave symmetry 
yielding a gap value Ao approaching closely onto the weak- 
coupling d-wave BCS relation Ao = 2.14IcbT c ; (2) At high 
doping levels, the residual SH term y(0) rises dramatically 
with doping, which suggests the existence of unpaired elec- 
trons possible in association with the normal metallic regions. 
These discoveries may carry out a common feature in cuprate 
superconductors and give important clues to the high-r c pair- 
ing mechanism. 
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